Observation of correlated optical and gamma emissions from GRB 081126 
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ABSTRACT 

We present an analysis of time-resolved optical emissions observed from the gamma-ray burst 
GRB 081126 during the prompt phase. The analysis employed time resolved photometry using opti- 
cal data obtained by the TAROT telescope, using BAT data from the Swift spacecraft and time resolved 
spectroscopy at high energies from the GBM instrument onboard the Fermi spacecraft. The optical emis- 
sion of GRB 081126 is found to be compatible with the second gamma emission pulse shifted by a 
positive time-lag of 8.4+3.9 sec. This is the first well resolved observation of a time lag between optical 
and gamma emissions during a gamma-ray burst. Our observations could potentially provide new con- 
straints on the fireball model for gamma ray burst early emissions. Furthermore, observations of time-lags 
between optical and gamma ray photons provides an exciting opportunity to constrain quantum gravity 
theories. 



Subject headings: gamma-ray : bursts 

1. Introduction 

Because of their extreme luminosity in y-rays, 
gamma ray bursts (GRBs) are a unique probe to high 
energy regimes where exotic physics is likely to mani- 
fest. A fraction of GRBs have been associated with the 
collapse of massive stars via the association of super- 



nova signature s observed with the fading GRB optica l 
afterglow e.g. faiorth et al.ll2003l: IStanek et al.ll2003b . 
The afterglow most Ukely originates from an exter- 
nal shock produced as the blast wave from the pro- 
genitor collides with the interstellar medium,causing 
it to slow down and lose ener gy. Fast moving tele- 



scopes linked to GCN notices (lBarthelmvlll998l) are 
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able to record the optical counterpart at the time when 
the prompt y-ray emission is still active. The first 
pos itive detec tion of such event was GRB 990123 
(lAkerlofetal.HT99'9i) . Some other successful detec- 



0.8)xlO"^ ergcm-2 



tions have been achieved so far (e.g. iRascusin et al. 
20081) . Two general results have been seen: Either a 



bright optical emission, uncorrelated to the gamma- 
ray light cur ve, occurred (for 5 to 20% of GRBs 



accor ding to iKlotz et al.1 l2009t) . or that a faint opti- 



cal emission is c orrelated with the ga mma-ray flares 
(GRB 050820A, I Vestrand et ai]|2006h . In the former 
case, these bright optical flash es are often inte rpreted 
as the reverse shock signature (iJin & Fanll2007h . 
Time lags between X-ray an d gamma-ray data are of- 
ten observed (e.g. iNorris et a l. 2000). However, this 
is r are between op t ical a nd y-rays. As an exam- 
ple, ^Tan g & Zhan 3 (12006') estimated the most prob- 
able time lags for the light curves of GRB 990123 
(5-7 sec) and GRB 0412 19A (1-5 sec). However, the 
optical data have poor time sampling, putting doubts 
on these results. Moreover, no la g was noticed for 
GRB 041219 bv lZheng etall Jiood) . No lag was re- 
ported for GRB 050 820A at a level of few seconds 
jvestrand et aDl2006h . 

In this letter, we present the measu rements of the op- 
tical emission observed by TAROT jKlotz et alj|2008l) 
during the prompt y-ray activity of GRB 081 126. We 
show evidence for a positive time lag between optical 
and y-ray light curves. 

2. GRB 081126 

GRB 081126 (Swift BAT trigger 335647, with 
To=26th Nov. 2008, 21:34:10 UT) light curve shows 
a small precursor starting at ~To-30 s, peaking at 
~ Tn-18 s, and r eturning almost to zero at T()-7 s 



(ISato et al.l l2008h . The burst features two peaks, 
the first one at ~To-i-1.5 s, reaching its maximum 
at ~To-h7 sec. The second one peaks at ~To+31.5 
sec. The duration of that burst is Tgo = 54 + 4 s (15-350 
keV). This event was a lso detected by Konus- Wind 
jGolenetskii et alJ l2008h and the Fermi GBM jBhatl 
2008h . 

The time-averaged spectrum of the first pulse from 
To to Tq-i-I 1 s is well fit by a Band function with Epggk 
= 192 + 74 keV, alpha = -0.3 + 0.4, and beta = -1.6 + 
0.1. The second pulse from about T()-h20 s to To+40 
s is also well fit by a Band function with Epcak -162 
± 11 keV, alpha = -0.3 + 0.5, and beta = -1.6 + 0.1. 
The fluence (8-1000 keV) in the two pulses are (2.7 



± 0.8)xlO"^ ergcm-2 and (1.9 
respectively. 

XRT observation reported by lMargutti et alJ (l2008h 
started 65.7 seconds after the BAT trigger, too late to 
gather X-ray information of the second peak. The 
XRT detected a characteristic afterglow emission of 
the burst. This afterglow was not detected by a quick 
visual inspection of ima ges taken by TARO T, started 



20.6 s after the burst jGendre et al.1 l2008h . How 



ever. lSkvarc & Mikuz ( 2008h reported the optical light 



curve of the afterglow in R band using the 60 cm of 
the Crni Vrh Observatory. Their observations start at 
T()+82s. They observe a slow rise in optical emission 
that peaks 200s after the trigger and then fade s. This 



optica l afterglow was also reported by Andreev et al. 



(12008D, using the Z-600 telescope of Mt. Terskol ob- 
servatory, 33 minutes after the burst, and by UVOT 
(Hollan d etal. 2008) at 2l''34"'03.59' -h48°42'38.3" 
(J2000.0). They report that the detection in the U fil- 
ter, combined with the lack of detections in the UV 
filters, is consistent with the afterglow having a red- 



shift of approximately 2.4 < z < 3.8 dHolland et al. 



2008 1). Unfortunately, no other photometric observa 
tions were performed to improve this estimation. From 
Konus-Wind data, we deduced a pse udo-redshift o f 
5.3 ± 1.8 using the method described in'A tteia ( 2003 ). 
The Galactic latitude of the afterglow position is -2.29° 
and the correspon ding extinction is E(B-V)=0.782 
mag. according to ISchlegel et al. I (fl998l) . Assuming 
.R=3.1, this gives Av=2.6 and A^=2.1 mag. 



3. TAROT data 

The first TAROT images were obtained at To -(-20. Is 
(duration 60 sec) with the tracking speed adapted to 
obtain a small trail of a few pixel length. This tech- 
nique is used in order to obtain temporal in formations 
during the exposure (e.g. lKlotz et al. I l2006h . The spa- 
tial sampling is 3.29 arcsec/pix and the FWHM of stars 
(in the perpendicular direction of the trail) is 2.05 pix- 
els. On the trailed image (see Figure[T]), the flux of the 
afterglow is affected by the proximity of NOMAD 1 
1387-0420537 (R= 18.1) but also by the end of the trail 
of NOMADl 1387-0420579 (R=15.48 hereafter A). 
This last star lies at 21 arcsec East and 7 arcsec South 
from the GRB position. As a matter of consequence, 
the trail of star A (which spreads over 30 arcsec) covers 
partially the beginning of the trail of the GRB (Fig.[T] 
top). 
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Knowing the position of the afterglow, we first sub- 
tracted the trail of the star A within the image. The 
star NOMADl 1387-0420302 (R=13.17, hereafter B) 
is far enough to other neighbor stars to be used as a 
trail template to model the star A. We then subtracted 
this model from the image (using a correct scaling fac- 
tor to take into account the difference of flux between 
the stars A and B). The result of the subtraction shows 
clearly the presence of a dim optical emission (Fig. [T] 
bottom). 



Successive images are 30 s long exposures tracked on 
iGendre et all (|2008|) pubUshed 



the diurnal motion, 
only upper limits using TAROT data because it was 
impossible to detect the optical counterpart so close 
to the star NOMADl 1387-0420537 without careful 
subtraction. The images taken later by TAROT were 
employed to perform this subtraction. The technique 
successfully revealed the optical afterglow. In Fig. |2] 
we display the initi al part of the TAROT lig ht curve. 
We add data from ISkvarc & Mikuz (l2008h showing 
that we can distinguish the early emission that oc- 
curred during the gamma activity and the afterglow 
that followed. A discussion of the afterglow emission 
process is beyond the scope of this paper and will be 
presented in Corsi et al. (2009, in preparation). 

4. Data analysis 

From the trailed image, a horizontal profile cor- 
responding to the predicted position of the afterglow 
gives directly the light curve. We measured the tempo- 
ral sampling of 6.5 sec/pixel using trails of bright stars. 
The light curve of the afterglow in the trail presents 
a flare within 3 pixels (Fig. O, with a probability of 
~ 10"^ to be spurious. The probability to observe a 
cosmic ray at that position is 3 .6 x 1 0"^ (estimated from 
dark fields of the same night). We thus conclude that 
this flare is real and produced by the burst itself. 

Such light curve profiles are affected by the Point 
spread function (PSF) of TAROT. In order to compare 
the optical and y-ray light curves, we need to con- 
volve the BAT signal by the TAROT PSF The PSF 
can be extracted as a vertical profile of a bright non 
saturated star (seen as the doted curve in the Figure|3j. 
We performed a symetrisation of the PSF shape to be 
compatible with the hypothesis that the PSF shows no 
direction effect. Once convolved with the TAROT PSF, 
the BAT signal corresponding to one peak is very sim- 



ilar to fliat of TAROT. 

We note at the start of the trail a bright single pixel 
that could be associated with the end of an optical 
flare. However, this event is not significant enough to 
be used in our analysis. Nonetheless, it could be an 
optical flare related to the first y-ray pulse. In the fol- 
lowing analysis, we will consider this as a possibility, 
and thus that the optical flare is correlated to the sec- 
ond y-ray pulse, without discarding the possibility that 
the optical flare is linked to the first y-ray pulse. 

A xl fit between the optical flare and y-ray pulse 
implies a temporal lag of +8.4 + 3.9 s (see Fig. |4]l at 
the 97% confidence (+38.4 + 3.9 s if the optical flare 
is related to the first y-ray pulse). This is strong evi- 
dence for a positive time-lag between the optical and 
high energy feature. We point out that the exposure 
time of TAROT images has a better accuracy than 0.1s 
because we use a GPS card triggered by the opening 
of the shutter, and is no t dependent on the com puter 
internal clock variations dLaas-Bourez et al ] E008l) . 

As the TAROT PSF is larger flian the BAT second 
pulse, we also studied the influence of the duration of 
the BAT pulse modelized by a Gaussian shape, let- 
ting free the width of the Gaussian within the fit. The 
best match of the modelized BAT pulse is a Gaussian 
spread by sigma=4.0sec. The xl fit gives the same 
lag as for the actual BAT pulse meaning that the pro- 
file shape of the pulse does not constrain the lag value. 
The fit remains compatible for Gaussians with sigma 
lower than 9 sec. This means the optical pulse is com- 
patible with a high energy pulse which could have a 
duration between to 9 sec. 

The flux of the optical peak observed by TAROT 
is 0.45 mJy. To be compared with the Fermi observa- 
tions, this value must be corrected for two effects: i) 
the spread of the flux due to the PSF profile, and ii) the 
large optical extinction in the R band. Correcting for 
all these effects, the optical flux is ~6 mJy at the peak. 
We used the Band model parameters obtained by the 
Fermi-GBM (Bhat 2008) to compute the optical flux 
expected from the high energy band. We derived an 
expected optical flux of 2.6x10 "*Jy, which is ~ 10"^ 
times the one observed. Taking account for the uncer- 
tainties in the Bhat's alpha parameter the extrapolated 
flux is always ~ 10"^ times the one observed. 

5. Discussion and conclusions 

The analysis of optical and gamma-ray light curves 
of GRB 081126 reveals: i) the width of the optical 
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peak is the same as the gamma-ray peaks, ii) the pro- 
file of the optical peak is consistent to the gamma- 
ray peaks after correcting for the different PSF, iii) 
the optical peak occurred 8.4+3.9 s (or +38.4 + 3.9 
s) later than the gamma peak. This is the first time- 
lag measured between optical and gamma light curves 
of a GRB. iv) the gamma-ray flux measured by GBM 
Fermi, extrapolated to optical energies is ~ 10"^ times 
smaller than the optical flux. 

These three results provide potentially new con- 
straints on the theory of prompt GRB emissions. For 
example, time lags between different energy photons 
are predicted by qua ntum gravity in the framewor k 



of string theory (e.g. lAmelino-Camelia et al.l Il998h . 



However, in such a case optical photons should arrive 
before gamma ones. As we observe the opposite, one 
can rule out this hypothesis for the GRB 081 126's op- 
tical lag. Gamma-ray photons comptonization on cold 
electrons could explain the profile of the optical flare. 
However, this cannot explain the positive lag observed. 

Within the internal shock framework, this temporal 
lag implies that optical photons were emitted after the 
y-ray ones. However, it is surprising that the flux in- 
creases so dramatically during this process. This is not 
well understood in the standard model for the inelas- 
tic internal shock and our results provide new tools for 
refining the standard model. 
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Eysseric, M. Merzougui, C. PoUas, P. Richaud and Y. 
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Fig. 1.— Field of GRB 081126. Top: TAROT im- 
age taken between 21s and 81s after the GRB trigger. 
The hour angle velocity was adapted to obtain stars 
as trails of ~9.2 pixel length during the 60s exposure. 
The theoretical position of the GRB trail is indicated 
by the white box. The star A covers partly the GRB 
trail, (see text). Bottom: After subtraction of star A 
using the model of star B, the trace of the optical emis- 
sion of the GRB appears in the box The image size is 
5 arcmin. North is up. East left. 
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Fig. 2.— Optical Hght curve of GRB 081126. 
TAR OT optical data a r e thic k bars and observations 
from ISkvarc & Mikua (120081) are thin bars (2 sigma 
level). There are data in the ranges 20-29 and 49- 
89 sec, but with no detection at the limit of R=18.0. 
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Fig. 3. — Light curves of GRB 081126 measured 
by BAT and TAROT. The dotted line labeled 'PSF- 
TAROT' stands for the spread of a star equivalent to 
an instantaneous flash of Os duration. This figure ap- 
pears in colors in the electronic version. 
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Fig. 4. — The convolution of the peak of BAT light 
curve (in blue) by the PSF-TAROT shifted by 8.4s (in 
black) compared to the TAROT data (in red). This fig- 
ure appears in colors in the electronic version. 
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